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Abstract. Monge-Ampere gravitation is a nonlinear modification of classical 
Newtonian gravitation, when the Monge-Ampere equation substitutes for the 
Poisson equation. We establish, through two applications of the large devi¬ 
ation principle, that the MA gravitation for a finite number of particles can 
be reduced, through a double application of the large deviation principle, to 
the simplest possible stochastic model: a collection of independent Brownian 
motions with vanishing noise. 


This paper is dedicated to Professor Tai-Ping Liu for his seventieth birthday. 

Introduction 

The purpose of this paper is twofold. We first want to make a short presentation 
of the Monge-Ampere gravitational (MAG) model, which has been introduced in 
m, in close connection with earlier works such as [n Ha El [1312] as well as with 
optimal transport theory (see Eill]). This model can be seen as a nonlinear mod¬ 
ification of the classical model of Newtonian gravitation, for which we use the fully 
nonlinear Monge-Ampere equation as a substitute for the linear Poisson equation 
to derive the gravitational potential from the density field. We also briefly com¬ 
pare the Monge-Ampere and the Newton gravitational models, emphasizing that in 
the MAG model there is an absolute control of the acceleration of the gravitating 
particles, no matter how concentrated they can be, in sharp contrast with Newto¬ 
nian gravitation, for which some particles may runaway to infinity in finite time. 
However, the main purpose of this paper is about deriving the MAG model from 
a completely elementary microscopic model in which a finite number of particles 
just move as independent Brownian trajectories without any interaction. In order 
to get the MAG model, we need two applications of the large deviation principle 
(LDP) [25l|21]. (Let us mention the interesting connection between large deviation 
principles, gradient flows and optimal transport theory recently made in 132] . 
which was certainly influential for us.) Through the first application of the LDP, 
we get a first order (in time) dynamical system based on the indistinguishability 
of the gravitating particles. (In a rather paradoxical way, the indistinguishability 
principle leads to a model where particles do interact.) The second application 
of the LDP enables us to lift this first order dynamical system to a second order 
one which is nothing but the discrete version of the MAG model. Let us emphasize 
that this derivation is purely formal and further investigations are clearly needed to 
get a complete and rigorous theory. The first section of this paper will be devoted 
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to a presentation of both the classical Newtonian and the Monge-Ampere models, 
with a suitable formulation of the MAG model relying on ” optimal transport” tools 
[iiiniisH]- The second section is devoted to the double application of the LDP in 
order to recover the discrete version of the MAG model. 


1. Monge-Ampere and Ne-wton gravitations 

1.1. Classical Newtonian gravitation. To describe classical Newton gravita¬ 
tion, let us attach to each gravitating particle a label a, which we suppose to belong 
to an abstract Borel measure space (A, A), and its position at time t, Xt{a) € 
(classically d = 3). For each time t, we denote by pt the image measure of A by Aj, 
defined on by 


(1) [ -/{x)pt{dx) = [ 7 (At(a))A(da) V 7 e C°(K‘^). 

For the sake of simplicity, we further assume pt{dx) to be periodic in x and 
of unit mean on the unit cube [ 0 , l]'^, so that we can see pt{dx) as a probability 
measure on the flat torus = K'^/Z'^. (This kind of assumption is common in 
computational Cosmology [26].) We are now ready to write the Newtonian model: 

( 2 ) =-VMMa)), = 


where the average density 1 has been subtracted out from the field pt so that the 
” gravitation potential” <pt is well defined from the Poisson equation as a zero mean 
Z'^—periodic function. (This is a perfectly meaningful assumption: there should be 
no resulting force out of a uniform pt .) If, at time t, pt is just a plain probability 
measure on then V(j)t, obtained through the Poisson equation A(j)t = pt — 1, 

is merely an function for p < (1 — on This makes very dubious the 

meaning of ODE 

d'^Xtia) 

—— — -X(pt{Xt{a)) 

in the general case. This is why Newtonian’s gravitation usually splits up into two, 
rather unrelated theories, according to the choice of its initial conditions: 

1) As A is a discrete set of N particles, then Newtonian gravitation boils down to 
an ”A-body problem”, set on the flat torus T'^, which is of form 


( 3 ) 


d^Xtja) 

dt"^ 


R{Xt{a)) + Y,G{Xt{a),Xt{b)), 

b^a 


where i? is a fixed function depending on the torus, and the singular kernel G{x, y), 
associated with the operator — VA“^, behaves as {y — x)/\y — x\'^, as x approaches 
y. This ’’bad” singularity makes the study of the ”A-body problem” very hard, 
in particular because of the possibility that particles may collide or runaway at 
infinity, in finite time. (See [18] as a recent reference.) 

2 ) As A corresponds to a ’’continuum” of particles, it is fruitful to introduce, at 
each time t, the nonnegative measure ft defined on the ’’phase space” T”* x Mf' by 

(4) / ^{x,v)Mdx,dv) = [ j{Xt{a), ^^^)X{da), V 7 G ^“(T'^ x R'^). 
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A straightforward calculation shows that / must satisfy, at least formally, the 
’’Vlasov-Poisson” system [35] 

(5) 9t/t + V, • (i;/t) - ■ (V0t/t) = 0, A(l>t = pt-l= f ft{dv) - 1. 

The meaning of this equation is unclear in the general case, because of the dubious 
product of V(/)t by /*. However, it is not hard to see that 0 < ft{dx, dv) < Cdxdv 
for some positive constant (7 is a consistent property for equation dSj). It turns out 
that this property is sufficient to justify the multiplication of ft by V(/)t at least 
when the initial condition /q vanishes for large values of v, say for |i;| > C for 
some other constant C. Then, the Newtonian model gets well defined and, at least 
for d < 3, global weak solutions can be shown to exist globally in time, as /o is 
given, whatever are constants C and C. In addition, such solutions are unique and 
smooth as long as /o has a smooth density with respect to the Lebesgue measure. 
(See [Ml [30l [8] for such results, stated on the whole space rather than on 
See also [35], as well as [22] for the closely related Euler-Poisson system.) 

1.2. Monge-Ampere gravitation. We now introduce the Monge-Ampere gravi¬ 
tation (MAG) model which differs from ([2]) just by the substitution of the Monge- 
Ampere equation for the Poisson equation: 

(6) = -VUMa)): det(I + D^t) = Pu 

where I denotes the dx d identity matrix and (ft is solution to the Monge-Ampere 
equation (in a suitable sense), which is periodic with zero mean and satisfies 
the (weak) ellipticity condition 1 + (j)t{x) > 0, in the sense of symmetric matrices, 

for every x. [Notice that (j)t is unique and Lipschitz continuous (resp. smooth) as 
soon as pt{dx) has an integrable (resp. smooth and positive) density with respect to 
the Lebesgue measure dx, see m for example.] We see that Newtonian gravitation 
can be formally recovered from the MAG model just by expanding the determinant 
about I and retaining only the linear part: 

det(I -I- D^4>t) ~ 1 -I- TraceD^^t = 1-1- 

Notice that, as d = 1 (which is a case of limited interest, describing gravitating 
“parallel pancakes”), the MAG model coincide with Newtonian gravitation. (As a 
consequence, our derivation of the MAG model from a double application of the 
large deviation principle, obtained in the second section of the present paper, is also 
valid for the Newtonian gravitational model in one space dimension.) However, let 
us emphasize that, to the best of our knowledge, MAG has never been considered 
by any physicist in dimension larger than one. It has been so far only considered 
by mathematicians (see m and [niniiiniia for closely related topics), mostly 
because of its close connection with optimal transport theory (as discussed in the 
next subsection). It is easy to describe, as we did for the Newtonian model, the 
MAG model through a ’’kinetic equation”, the so-called ’’Vlasov-Monge-Ampere” 
(VMA) system [TSII5T] : 


(7) dtft + y.-{vft)-yv-{^^tft) = 0, det{l +D^<Pt)=Pt 


ft{dv). 
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As for the Vlasov-Poisson system the existence of global weak solutions can be 
shown as soon as 

0 < fo{dx, dv) < Cdxdv, fo{T‘^ x {u e \v\ < C'}) = 0, 

for some positive constants C,C' . Existence of a unique smooth solution, but only 
for a short time interval, has been proven by Loeper m provided fo{dx,dv) (resp. 
po{dx)) has a smooth density with respect to dxdv (resp. dx) and is uniformly 
compactly supported in u € (resp. is strictly positive). These are still limited 
results, due to the full non-linearity of the Monge-Ampere equation that leads to 
analytic difficulties, in particular when seeking for smooth solutions. However, the 
Monge-Ampere equation enjoys remarkable properties, closely related to the theory 
of optimal transportation [iniiM]- This is why we are going to introduce a related 
formulation of the MAG model, with interesting geometric features. 

1.3. The MAG model written in optimal transportation terms. For the 

description of the MAG model in optimal transport terms, it is convenient to discuss 
the model in a (slightly) different and more abstract framework. We consider a 
metric measure space (A, A), made of a compact subset A of equipped with a 
Borel probability measure A. Two typical examples are, on one side, the unit cube 
with the Lebesgue measure, and, on the other side, any set of N points equipped 
with the (normalized) counting measure. In the first case, we will speak of the 
’’continuous” case, while, in the second case, we will speak of the ’’discrete” case. 
(Of course many others situations could be also considered, in particular the flat 
torus T'’* as we did in the previous subsections, but we will focus on these two 
cases only.) We introduce the separable Hilbert space H of all A—square-integrable 
maps from A to K'^, H = L^(A, A;K'^), with norm and inner product respectively 
denoted by || • || and ((•,•))■ (Notice that, when A is a finite set of N points in 
then H ^ is of finite dimension.) We crucially consider the subset S of all 
measure-preserving maps s of A: 

( 8 ) S = {s€H, [ 7 (s(a))A(da) = [ 7 (a)A(da), V 7 e C'°(R'^)}. 

Ja Ja 

(Observe, in the discrete case, when A is made of N distinct points A(a) G R'^, for 
a = 1, - ■ ■ , N, S just reads S = {(A(cr(l)), • • • , A{a{N)) € H, a € Sn} with N\ 
elements, where Sjy denotes the group of all permutations of the first N integers.) 

According to Edelstein’s theorem [ 23111 ], given a separable Hilbert space H and 
a closed bounded subset S, there is, in the sense of Baire, a generic set (i.e. con¬ 
taining a countable intersection of dense open subsets of H) of points X for which 
there exists a unique closest point 7r(A) on S. In addition 7r(A) is nothing but 
the gradient, at point X, of the Lipschitz convex function H defined on the Hilbert 
space H by 

(9) n(X)=sup((X,s))-K, ^(x)=Vn(X)=Arginf 

ses 4 ses z 

In our particular case, we can say much more in the continuous case, when A is the 
unit cube with A as the Lebesgue measure, thanks to ’’optimal transport theory” 

Theorem 1.1. /^[lOj ) Let X G H be a non degenerate map, in the sense that the 
image measure p of X by X is absolutely continuous with respect to the Lebesgue 
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measure on Then, X has a unique closest point tt{X) on S. We also have 

(10) tt{X)=ToX, 

where T is the unique map (in the p a.e. sense) y G —>■ T(y) ^ A such that: 

i) there is a convex Lipschitz function ■i/i : ^ R, p—a.e. differentiable, with: 

(11) T{y)=Xip{y), p-a.e.y 

ii) X is the image of p by T: 

(12) [ ^{T{y))p{dy)= [ i{a)X{da), V 7 G C°(R‘'). 

JRl JA 

As a matter of fact, equations (I11I12I) form a generalized formulation of the 
Monge-Ampere problem on R'^ 

(13) p = det(iD^i/;), (Vi/')(support(p)) = A, > 0, 

with a unique solution Xif (in the p—a.e. sense). (See [9l[Tnil38] for more details.) 
This generalized formulation of the Monge-Ampere equation allows us to write the 
MAG model ([HI) in a much more geometric way. Indeed, m just reads 

= Xt{a) - V-iftiXtia)), det{D'^tjjt) = pt, 

after setting ipt{x) = \x(^/2-\-(j)t{x), that we complete with the weak ellipticity con¬ 
dition > 0 and the range condition V'!/'t(support(pt)) = A (which substitutes 
for the Z'^-periodicity condition we used in writing (O). Assume that, at time t, 
Xt is non degenerate (or, in other words, pt is absolutely continuous with respect 
to the Lebesgue measure on R'^). Then, using Theorem ll.il we may write 

Vfjt o Xt = Tr{Xt) = XU{Xt) 

and finally obtain 

(14) ^^=Xt-7roXt = Xt-XUoXt, 

where II is the Lipschitz convex function defined by ®, TT its gradient and S is the 
set of measure preserving maps defined by (jS]). In the rest of this paper, we will 
retain (I8I9I14P as our definition of the MAG model. 

1.4. The discrete Monge-Ampere gravitational model. The geometric for¬ 
mulation (I8I9I14I) of the MAG model is very convenient to get its discrete version, 
as A is just a subset of N points A(a) in R'^, for a = 1 , • • • , A, equipped with the 
counting measure, in which case, we will speak of the ’’discrete MAG model with N 
particles”. Indeed, at the discrete level, a time-dependent map Xt can be seen as 
a sequence of N ’’particles”, with positions Xt{a), moving in R'^, for a = 1, • • • ,N. 
(Notice that Xt gets degenerate at time t in case of collisions.) In this discrete 
setting, the MAG model (1819114^ reads as the finite dimensional dynamical system: 

72 V ( \ 

= ^t(a) - ^(iTt(a)), crt = Arg inf ^ |At(a) - A(CT(a))p, 

cre5jv 

a—\ 

where a G Sjs; denotes the set of all permutations of the first N integers. So we 
see that, at the discrete level, the MAG looks both very simple and very different 
from the classical Newtonian N body problem ([3])! (Similar systems have been 
previously studied in mm) 
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The MAG model, as defined by (I8I9I14I1 . can be seen as a dynamical system in 
a Hilbert space H with a force term J^{x) = x — Vn(a;) with a trivial linear part 
and the gradient of a Lipschitz convex function H. To the best of our knowledge, 
there is no theory for such an ODE in infinite dimension. Indeed, we are very far 
from the standard Cauchy-Lipschitz setting. This is a challenging open problem. 
(See the related theory developed by Ambrosio and Gangbo for some infinite di¬ 
mensional hamiltonian systems [5]. As a matter of fact, their main example is very 
similar to the MAG model, written in a different way. See also [13 HZ].) However, 
for the discrete MAG model with N particles, H ^ is of finite dimension. 
Thus, equation (fT^ . which also reads m, can be neatly solved globally in time 
in the sense of Bouchut and Ambrosio mm, for every initial condition (AIq, 
except on a negligible subset of the ’’phase space” H x H, oi zero 2Nd dimensional 
Lebesgue measure. (A more accurate statement can be found in [T].) Indeed, the 
force term T(x) = x — Vn(a:) in the right-hand side of (fl3 is a smooth pertur¬ 
bation of a ’’bounded variation” function, since H is a Lipschitz convex function 
(which implies that D^Il can be seen as a bounded nonnegative measure valued in 
the convex cone of nonnegative symmetric matrices, and, therefore, that VH is of 
bounded variation). Notice that the exceptional set of bad initial conditions in the 
phase space is not empty, as it follows clearly from formulation (nsi), where we see 
that the evolution of the particles becomes ambiguous as different particles depart 
from the same position with exactly the same velocity. 

1.5. Monge-Ampere versus Newtonian gravitational models. It is now in¬ 
teresting to compare the Monge-Ampere and the Newton gravitational models. We 
first observe that, according to the MAG model, particles may never runaway to 
infinity in finite time. Indeed, from the optimal transportation formulation (I8I9I14L 
we immediately get 

J2 V 

I^-WI < i?=sup H 
oeA 

and deduce that 

\^t\ + \~^\ — (l^ol + I I + 1)C cosht, Vt G R, 

where C depends only on R. Of course, there is nothing similar with Newtonian 
gravitation. Indeed, the Poisson equation in ([2]) does not behave well when particles 
concentrate. For instance, if some particles concentrate as a delta measure at some 
point y at time t, then V(j)t{x) has a singularity as bad as {x — y)/\x — y\‘^. 

In this way, the MAG model has a lot of similarity with the Born-Infeld (BI) theory 
of the electromagnetic field [3 Emmie], in which any electrostatic field is bounded 
by a universal constant. (Notice that the BI model, which goes back to 1934, is no 
longer for use in Electrodynamics, but has enjoyed a remarkable revival in String 
Theory since the 1990s [ 33 ] •) 

In addition, the MAG model enjoys good properties of approximations by finite 
sums of Dirac measures (see [20] for closely related results). In sharp contrast, such 
discrete approximations have never been justified, to the best of our knowledge, 
in the case of 3D Newtonian gravitation, because of the bad singularities of the 
Green function for the Poisson equation [28] ■ As a matter of fact, the treatment of 
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point particles in classical Electrodynamics remains an outstanding open problem 
in both Theoretical and Mathematical Physics (see [24l [37] for instance). 


2. Formal derivation of the discrete Monge-Ampere gravitation model 
from a donble application of the large deviation principle 

In this second section, we provide a formal derivation of the discrete MAG 
model ([T5|1 from a double application of the large deviation principle [25l [21]. 

2.1. A basic model of N independent Brownian particles. We fix a positive 
integer N and a finite set A of A distinct points A{a) G for a = 1, • • • , A. 
Using the notations of subsection 11.31 we set A = and 

(16) S={(A(ct( 1)),-- - ,A(a(N))eII, a € Sjv}, 

where Sjy denotes the set of all possible permutations of the first A integers. 
Given ctq G Sn, we consider the motion of A ’’particles” that move in the Euclidean 
space according to 

(17) X^{a) = So{a) +eBt{a), So{a) = A{ao{a)), Va = l,---,A. 

where {t G R+ —S' Bt{a) G R‘^)a=i, -.v denote A independent realizations of the 
’’standard” (i.e. normalized) Brownian motion in R'’*. 

Fixing U > 0 and 

F*(a) eR'^, Va= I,-- - ,A, 

it is easy to compute the probability that, at time t*, each particle occupies the 
position given by Y*-. 

Prob(Af. (a)s»y*(a), Va = 1, • • • , A) 

a—1 

or, in other words 


Prob(Af.^y*) exp(^^-^^)(27r£G)-^'^/2. 

Here we have denoted by | • | and 11 -11 the Euclidean distance on respectively R*^ and 
H = (R‘^)'^. We have also used symbol Ri just to make notations lighter. [What 
we precisely mean is: for any Borel B subset of (R*^)-^, the probability that Af. 
belongs toY*+B is given by 

/ exp(^t^)(2^er)-""/2dy, 

JY^y+B 


but we hope that our simplified notation is acceptable. Indeed, we will use it again, 
without further notice.] 
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2.2. First application of the large deviation principle. We now want to com¬ 
pute the probability that Xft and Y* coincide, up to a permutation, property that 
we denote by (X^,) fa {¥*). We find 

perm 

(18) Prob[(Xf.) ^ (F*)] « ^ 

perm Jy ! ^^ ZSt* 

(tGSn 


When the level of noise e goes to zero, we immediately get 


— lim elogProb[(X|,) (F*)]ss inf 

£—)-0 perm <tGSn 


\\Y* o cr - SqIP 
2t* 


which is a rather trivial illustration of the Laplace method and the large deviation 
principe (LDP) [^. There is a more sophisticated aspect of the LDP: as the 
level of noise goes to zero, the Brownian trajectories of the particles conditioned 
by (Xft) fa (Y*) behave more and more as constant speed minimizing geodesic 

perm 

curves. More precisely, for all t S [0,t*], 


t IIF 

(19) X^ Xt = So + —{Y* oa* - So), a* = Arg inf J-— 

L Zii 


o cr — Sol 


which implies 




F* O CT* - So Xt* - So Xt- So 


vt e]o,r 


t* t* t 

Now, we are going to translate this large deviation principle into a self-consistent 
dynamical system for the particles. Let us first denote by 7r(F) (as in subsection 
fT31) the unique closest point on 5* of a "generic” point Y G H = 


N 


( 21 ) 


7r(F) = Arg inf 


||F-s|p 


= Vn(F), n(F) = sup ((F, s)) - ^ 

seS 2 


(where ((•, •)) denotes the inner product attached to || • || on H). Then we state: 


Proposition 2.1. Equation i2CKl . that we have derived from (through a large 
deviation principle), can be written as a self-consistent ordinary differential equation 
for Xt at least for t > 0, 

(22) t^ =Xt- Tr{Xt). 

Proof. To get this result, our simple but crucial observation is that, along the 
geodesic curve defined by (HU), namely 

At = So-f ;^(F* o cr* - So), cr*=Arg inf ||F* o cr - so||. 

So is the closest point in S not only of the end-point Xt* = Y* o a* but also of all 
points Xt, Vt G [0,t*]. [Then, we can write sq = 7 r(At) in equation (l20l) which 
immediately leads to (1^^ and completes the proof of our Proposition.] 

Although this property is geometrically quite obvious, let us provide a comprehen¬ 
sive proof for the sake of completeness. 

By definition m of the closest point operator tt, it is enough to show that 
K = ||At - So o crjp - ||A( - SolP > 0 
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for each permutation a G Sn and each t G By definition (fT^ . 

= ll'So + “ ®o) - So o (t||^ - W^iX* - so)|p 

where we set, still according to (HU), 

(23) X*=Xt.=Y*oa\ cr* = Arg inf ||r* OCT-Soil- 

(t^Sn 

Expanding the squares, we obtain 

2t 

K = ||so - So O cr|p + —{{X* - So, So - So O aj). 

By definition of X* =Y* o a* and a*, we get from (1^51) 

||X* - Soil = ||y* OCT* - Soil < ||y* 0(T*0CT-1 - Soil = ||A:* -SoOct||. 

Thus ((X*, So — So o a)) > 0. So, we deduce 

2t t 

K> ||so - So o cr|p - —((so. So - So o cr)) = 2(1 - —)[||so||^ - ((so, So O cr))] > 0 

(using 11 So 11 = 11 So o (T11 and the Cauchy-Schwarz inequality), which is the desired 
inequality and completes the proof of Proposition 12.II 

Equation (1^^ is clearly singular at time t = 0. We can lift this singularity with an 
exponential rescaling of time t = exp(0), 0 G M, and hnally obtain: 

(24) ^=Xe- TT{Xe) = Xg - Vn(Xe). 

dO 

So we have obtained a first order (in time) dynamical system, of ” gradient type”. 

Miscellaneous remarks, i) Since tt is the gradient of a Lipschitz convex function 
(namely 11 defined by ®), equation (IMl) is uniquely solvable in the framework of 
’’maximal monotone operator theory” m- As already discussed in m, this is a 
way to introduce a dissipative mechanism in the motion of particles, such as sticky 
collisions when particles stick to each other while conserving their momentum (but 
not their kinetic energy which decreases). At the level of the present paper, we do 
not want to enter such considerations and leave (1241) just as a formal equation. 

ii) Notice that (1^^ is a gradient flow in the variable X, valued in the Hilbert space 
H, which is the counterpart of the gradient flow in the variable Law(Ai), in the 
so-called ’’Wasserstein space” of half the negative squared ’’Wasserstein distance”, 
as discussed in [3]. For us, it is important to keep a formulation in terms oi X G H 
and not in terms of Law(X), because, in the sequel of our discussion, it is crucial 
to restore the ’’individuality” of the gravitating particles during their motion, in 
order to get a second order dynamical system for them. 

iii) Quite remarkably, as explained in m , equation (l22l) is nothing but the Zeldovich 
model used in Cosmology [311 [Ml Hi [H as an approximation of semi-Newtonian 
gravitation in an Einstein-de Sitter space! 

iv) Let us finally provide a possible interpretation of equation (1^ , vaguely related 
to the so-called ’’anthropic principle”: What is observed at out “present” time 
t* is just a random output of the independent Brownian trajectories of a large 
number N of indistinguishable particles initially located on the set A. As the 
noise vanishes, the motion of these particles, conditioned by what we can observe 
now, just looks driven by the deterministic law (12411 . Of course, this is a highly 
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speculative and questionable interpretation coming from a mathematician and not 
from an authorized physicist. 

2.3. Second application of the large deviation principle. Equation (1241) is 
a first order dynamical system, a so-called ’’gradient flow” since ((24ll can also be 
written 

dXs = V^{Xg)de, 

where 4) is half of the squared distance function to S: 

(25) $(y) = ffl_n(y) = ig^l^, 

where If is defined by ([9]). Describing a gravitational theory by a gradient flow does 
not sound reasonable. We would rather like to get a second order, conservative, 
dynamical system. Here again, the large deviation principle turns out to be useful 
to get such a second order system out of (l24ll . Fixing 77 > 0, we introduce the 
’’noisy” version of (1^51) defined by: 

(26) dX^ = X^{X^)de + T] dBg, 

where Bg (again) denotes a Brownian process in . Given two points Id and 

Yi in , we expect from the LDP (or, more precisely, the Freidlin-Wentzell 

theorem [251 [2T]) that, as the level of noise rj goes to zero, 

-jylogProb[X^^ssYo andX^^RiYi] -^^0 A{9o,Yo,9i,Yi) 

(27) A{9o,Yo,9i,Yi) =M{ r -X^Xg)\\^d9, Xg^=Yo, Xg,=Y,}. 

Jbq ^ 

In addition, as r] goes to zero, 

XI ^Xg, V0e[0o,0i], 

where 

X = Arginf{^''I||^-V4>(Xe)||2d0, = Yq, = Yi}. 

Strictly speaking, this is correct when V$ is Lipschitz continuous, which is not true 
in our case (where V<i> is not even continuous), without further assumptions on the 
data. However, from the pure modeling viewpoint, it is very tempting to find the 
second order dynamical system linked to the least action principle (I27p . Since $, 
as defined by (EH), is half of a squared distance function, we have 

(28) i||V4>(Y)||2 = ci>(Y), 

for every Y € H \ JV, where Af is the set on which $ is not differentiable, which is 
a negligible subset of H, both in the Lebesgue almost everywhere sense and in the 
Baire category sense. Thus, at least for each curve X that stays away from Af for 
Lebesgue almost every 9 G [ 00 ,0i], we have 

C 7 "^ - 

= + l|V4>(A,)ind0 - 1 '^ y^Xg) ■ ^d9 
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= + HXe))d9 - 

Therefore, the action principle (12711 is equivalent to 

(29) Ai9o,Yo,9uYi)=ini{j\^\\^\\^ + ^Xe))de, Xe,=Yo, Xs,=Yi}. 
From this equivalent LAP, we find as optimality equation 

(30) ^ = V<i>(X,) = Ae - iriXe), 

which is just the second order version of the gradient flow equation (1241) . More 
explicitly, we have obtained 

(31) ^^j^ = Xo{a)-A{ao{a)), 

where 

N 

ere = Arg inf V |Ae(a) - A(cr(a))p, 

(tGSn 

a—1 

which is nothing but the discrete version (HID of the Monge-Ampere model of grav¬ 
itation discussed in the first section. 

So we have achieved, at least at a formal level, the main goal of our paper which 
was, through a double application of the large deviation principle, the derivation of 
the (discrete) Monge-Ampere gravitational model from one of the simplest think¬ 
able model of particles: N independent Brownian trajectories, with an intriguing 
interplay between their indistinguishability (for the first application of the LDP) 
and their individuality (for the second application of the LDP). 
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